Saturation of electrical resistivity in metals at large temperatures 
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We present a microscopic model for systems showing resistivity saturation. An essentially exact 
quantum Monte-Carlo calculation demonstrates that the model describes saturation. We give a sim- 
ple explanation for saturation, using charge conservation and considering the limit where thermally 
excited phonons have destroyed the periodicity. Crucial model features are phonons coupling to the 
hopping matrix elements and a unit cell with several atoms. We demonstrate the difference to a 
model of alkali-doped Ceo with coupling to the level positions, for which there is no saturation. 



In a metal, the electrical resistivity, p, grows with the 
temperature T due to the increased scattering of the elec- 
trons by phonons. Typically, p{T) ~ T for large T. For 
some metals with a very large p, however, the resistivity 
saturates |^-|^, i.e., it grows very slowly with T for large 
T. The resistivity is often described in a semiclassical 
(Boltzmann) picture, where an electron, on the average, 
travels a mean free path I before it is scattered. The re- 
sistivity is inversely proportional to /. Typically, I 3> d, 
where d is the atomic separation. For systems with re- 
sistivity saturation, however, I becomes comparable to d. 
Work in the 1970's suggested that resistivity saturation 
occurs universally when I ^ d, the loffe-Regel condition 
Q , providing an upper limit to the large T resistivity of 
metals. Later work has, however, found exceptions, such 
as alkali-doped fuUerenes |^,^ . 

Intuitively, resistivity saturation seems natural. One 
might expect that at worst, an electron could be scattered 
at each atom, leading to / d. Such a semiclassical pic- 
ture, however, breaks down when ^ ~ d 0, and it is con- 
tradicted by the lack of saturation for fuUerenes. Several 
theories of the saturation have been presented, usually 
based on generalisations of the semiclassical Boltzmann 
theory, but none has been generally accepted [pj-p^. 

The Bloch-Boltzmann theory starts from a periodic 
system and treats the scattering mechanisms as small 
perturbations. Here we consider the opposite limit, 
where thermally excited phonons have removed all effects 
of periodicity. In this limit, charge conservation natu- 
rally leads to saturation for systems with several atoms 
per unit cell and strong electron-phonon coupling, in par- 
ticular for systems where the phonons couple to the hop- 
ping matrix elements. We show that this happens to 
occur when I ~ d. This does not happen, however, for a 
model of alkali-doped Ceo, where the phonons couple to 
the level positions. We first use an essentially exact quan- 
tum Monte-Carlo (QMC) calculation to demonstrate sat- 
uration in our model. We then introduce a method where 
the phonons are treated semiclassically and the electrons 
quantum-mechanically, justifying the method by compar- 
ing with the QMC results. This method is sufficiently 
simple to allow for an interpretation of the results. 

Saturation is clearly observed for, e.g., A15 com- 
pounds, such as NbaSb [Q, while Nb metal shows weak 



saturation at large T ||Tl|. We study a model Nbg of 
NbsSb, where the Nb atoms have the same positions as 
in NbaSb, but the Sb atoms are neglected This is 

compared with a model of Nb. We consider a cluster of 
N atoms, placed on A15 (Nb3) or bcc (Nb) lattices with 
the lattice parameters 5.17 A (Nbg) or 3.28 A (Nb). 

We study the scattering of the electrons from phonons. 
Each atom is assumed to have vibrations in the three co- 
ordinate directions, described by Einstein phonons. The 
phonon energy, cuph = 14 meV, is obtained from an aver- 
age over the phonons of Nb [|l3|. For each Nb atom we 
include the five-fold degenerate (n = 5) c?-orbital. The 
hopping matrix elements between the orbitals are ob- 
tained from Harrison ]l^ , using a power dependence on 
the atomic distance d (~ l/c?™). We use m = 3.6 more 
appropriate for Nb ||l5| than to = 5 used by Harrison. 
To avoid divergencies for very small d, 1/d™ is replaced 
by l/(d™ + a™), with a — 2 k. The atomic vibrations 
modulate the hopping matrix elements, both due to the 
changes of the atomic distances and the relative move- 
ments of the orbital lobes. We neglect the influence of the 
vibrations on the level energies as well as the Coulomb 
interaction between the electrons. To obtain the resistiv- 
ity, we calculate the current-current correlation function. 
The current operator j is obtained by using charge and 
current conservation. Thus the matrix elements between 
orbitals ly and /i on sites R'' and are 

r''=ie(R''-R'')t,.^, (1) 



where t^^ are the corresponding hopping matrix ele- 
ments. 

This model can be solved essentially exactly by using 
a determinantal quantum Monte-Carlo (QMC) method 
Jl6t , treating the phonons quantum mechanically. For 
the models studied here, the QMC method has no "sign- 
problem" . The calculated correlation function (for imag- 
inary time) therefore has only (small) statistical errors. 
We use a maximum entropy method [ p^ to obtain the op- 
tical conductivity (t(lu) on the real frequency axis. The 
resistivity p is then l/a{0). 
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FIG. 1. Resistivity /9(T) as function of temperature T. The 
main figure shows resuhs for Nb3 (A'^ — 648 atoms) and the 
right insert for Nb (A*' = 640 atoms). In both cases the 
phonons are treated semiclassically. The small T (Eq. (Q)) 
(chain curves) and the large T (Eq. ^) (broken curves) lim- 
its of p(r) are also shown. The left insert compares the semi- 
classical (broken (A^ = 36) and full (A*' = 648) curves) and 
quantum Monte-Carlo (circles, N = 36) calculations for NbJ^. 

The left insert of Fig. |^ shows the QMC result (circles) 
for p{T) of the model. We are particularly interested 
in the large T behaviour, which is also the limit where 
the QMC calculation can be performed with a reasonable 
effort. The large T result extrapolates to a substantial 
nonzero value. However, since the resistivity of the Nb3 
model must go to zero for T = 0, the QMC calculation 
clearly shows that the model leads to a drastic reduction 
of the slope of the p{T) curve for large T, usually referred 
to as resistivity saturation. 

To analyse the results, we treat the phonons 
(semi)classically. The atomic displacements due to the 
phonons are chosen randomly according to a Gaussian 
distribution, determined from the average number of 
phonons at that temperature. For given "frozen" dis- 
placements, we calculate the hopping and current matrix 
elements. The eigenvalues Ei and eigenvectors \i) of the 
resulting Hamiltonian are calculated. For an isotropic 
system, the optical conductivity is then given by 

a{cu) ^ i J2\^^\UJ)\^if^ - fj)S{hcj ~ e, + e.), (2) 

ij 

where fi is the Fermi function for the state i. 

The left insert of Fig. Q shows that this approach (bro- 
ken line) agrees quite well with the QMC calculation for 
large T, the temperature range we are interested in, and 
we expect the semiclassical calculation to remain accu- 
rate for T ojp/i (=0.014 eV). We can therefore interpret 
the results by analysing the simpler semiclassical calcu- 
lation. We observe that this model differs qualitatively 
from the (Ziman solution ||l^) of the Boltzmann equa- 
tion, in which p(T) ~ T for large T. In the Boltzmann 
equation the electrons are treated semiclassically, while 



here they are treated quantum-mechanically. 

Fig. |l| shows that there is a very clear saturation for 
NbJ, while for Nb (right insert) there is only weak satu- 
ration at large T, in agreement with experiment |p],pT|. 
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FIG. 2. Optical conductivity as function of frequency uj for 
the (a) A15 and (b) fuUerene models in the semiclassical cal- 
culation. The frequency has been scaled by the T = band 
width W . (a) also shows (broken curve) the result of approx- 
imating all current matrix elements by their average. 

We now discuss this saturation. Fig. ||a shows a{oj) 
for Nb3 in the semiclassical theory. For small T, there 
is a narrow Drude peak at w = 0, which is smeared out 
for large T. Then a{io = 0) drops correspondingly. For 
u! > W, (j{lu) remains zero (apart from a slight broaden- 
ing introduced in the calculation), where W is the band 
width, since there are no excitations for u > W . 

The Drude peak is due to intraband transitions be- 
tween states with similar values of the wave vector k. As 
T is raised and the vibrations of the atoms are increased, 
the states loose their well-defined k and band index la- 
bels. To illustrate this, we decompose a state \i) at a 
finite T in the T — states with given k-vectors. The 
corresponding weights are labelled c(ik). We define A as 
the average of A{i) over all states i, where 

A(*)=nk^c(ik)2, (3) 

k 

and rik = 256 is the number of allowed k-vectors for the 
systems studied. If each state has rik^/m k-components 
with equal weights, A ~ m. If periodicity is completely 
lost, A = 1, and if each state has only one k-component, 
A = nk(= 256) A is shown in Fig. |. It illustrates 
how the effects of periodicity are lost very quickly for Nbg 
on a temperature scale of just a few hundred K. 

For a large T, the (q 0) current operator then cou- 
ples all states to each others. We now make the assump- 
tion that at large T the coupling between all states is 
equally strong This is the opposite limit to the 

Boltzmann treatment, where k is assumed to be a good 
quantum number. We replace a current matrix element 
in Eq. (||) by its average over all transitions. 
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FIG. 3. The measure A (Eq. (Bh) of the amount of momen- 
tum conservation as a function of T for Nbs (full line) and Nb 
(broken line). The horizontal line represents complete loss of 
momentum conservation. 

Using charge and current conservation (Eq. (|l])), we ob- 
tain 



3Nn- 
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(4) 



where n is the orbital degeneracy, (e^) is the second mo- 
ment of the density of states N{e) per atom and d is 
a typical atomic separation. We assume some generic 
shape of N{e). Since we know its second moment, we 
can then relate N{e) to the j^^. Using this relation, we 
can calculate cr(w), shown by the broken curve in Fig. 

A comparision with the full calculation shows that 
this approximation becomes rather accurate already at a 
moderate T. We find that cr(0) ^ (P/^l, where is the 
volume per atom. Using the T = nearest neighbour 
distance for d, we obtain 



n 



(5) 



A depends on the precise assumptions of the model, but 
it is ~ 0.2 mficm if d is expressed in A. Assuming a semi- 
elliptical density of states ~ ^ {W/2)'^ — and band fill- 
ing 0.4, we find A = 0.27 and 0.15 milcm for the A15 and 
bcc lattices, respectively. This is shown by the broken 
horizontal lines in Fig. |l|. Given the simple assumptions, 
the agreement with the full semiclassical calculations at 
large T is surprisingly good. The corresponding apparent 
mean-free path is 



I 



i^'^d. 



(6) 



With the assumptions above, we obtain c — 0.5 and 
c — 0.6 for the A15 and bcc lattices, respectively. Sat- 
uration therefore happens roughly when the loffe-Regel 
condition is satisfied, as might have been expected on 
dimensional grounds. 

The derivation of Eq. (^ uses a quantum-mechanical 
treatment of the electrons. It explains why metals with 



a large resistivity usually show saturation, and why it 
happens when I ^ d. 

The p(T) in Eq. (|) is independent of T, while Fig. 1 
shows a weak T dependence even after saturation. This 
is partly due to the T-dependence of the Fermi-functions 
in Eq. (^, which were approximated by 0- functions in 
the derivation above. There is also a T dependence due 
to changes of the shape of the density of states as T is in- 
creased, which was neglected above by using a "generic" 
density of states. These effects are only important when 
T is an appreciable fraction of the band width (for real- 
istic values of the electron-phonon coupling) . 

Eqs. (H, ^ are obtained by assuming that the Drude 
peak has been smeared out and that <t{uj) is spread out 
over the whole band width. Both the saturation resistiv- 
ity and the corresponding mean free path are indepen- 
dent of the band width. A scaling of by a factor a > 1 
increases W, and o'(ijj) extends over a larger energy range. 
At the same time, however, j^^ is increased (Eq. (Q)) in 
such a way that cr(0) in Eq. (||) is unchanged. Charge 
and current conservation therefore plays a crucial role for 
our results (|,|) . 

We next consider a small r(> utph). Then 



piT) = Stt- 



, XTkB 



(7) 



where A is the electron-phonon coupling constant and fc^ 
is the Boltzmann constant. Qpi is the plasma frequency 



Ef). (8) 



where £„k is the energy of a state with the band index n 
and the wave vector k and Ep is the Fermi energy, flpi 
depends on the average Fermi velocity. 

The straight line given by Eq. (^ (chain lines in 
Fig. 0), agrees well with the semiclassical calculations 
for small T. Typically, this line rises so slowly that it 
intersects the horizontal line of Eq. (^) well above the 
melting point for the metal of interest. Then no satu- 
ration is found in the accessible temperature range. If, 
however, flpi is small and A is large, the two lines cross 
below the melting temperature. Saturation is then ob- 
served. An important difference between Nbg and Nb is 
that Nb3 has many rather flat bands, due to the large 
unit cell |^ and many forbidden crossings. The result- 
ing small electron velocities, lead to a small ilpi . We find 
^Ipi — 3.6 and 8.2 eV for Nb3 and Nb, respectively, which 
makes the slope of the line in Eq. (|^) about a factor of 
five larger for Nbg than for Nb. We find similar A's for 
Nb;^ (A ^ 1.0) and Nb (A = 0.9). More accurate esti- 
mates give npi = 3.4 eV (NbgSn) Q and 9.5 eV (Nb) 
H and A = 1.7 (NbgSn) [§ and 1.1 (Nb) 

It is interesting to replace the d-orbitals in our Nbg 
model by s-orbitals. As before, the resistivity shows sat- 
uration, but the saturation is a less pronounced than in 
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Fig. 1. The saturation resistivity (Eq. (|^)) is larger due 
to the smaher degeneracy (n = 1). 

A very different behaviour is found in the alkah-doped 
fuUerenes, (AaCgo (A= K, Rb)), where the apparent 
mean free path becomes much shorter than the sepa- 
rations of the Ceo molecules |^,||. This behaviour was 
related to the fact that the important (intramolecular) 
phonons primarily couple to the Ceo level energies, in- 
stead of the hopping matrix elements | p3| , p^ . Already 
at a moderate T, the resulting fluctuations in the level 
energies become comparable to the T = width of the 
narrow tiu band, which conducts the current. This leads 
to a broadening of the tiu band and of cr(w) beyond the 
T = band width. In contrast to the case of a scaling of 
the hopping parameters, discussed below Eq. (^, this is, 
however, not accompanied by an increase of the current 
operator. Thus a{Q) is reduced, explaining the lack of 
saturation. This is illustrated in Fig. ^d. 

It is interesting to compare the cases when the phonons 
couple to the level energies (LE coupling) and to the hop- 
ping integrals (HI coupling). We have studied the case 
of HI coupling in a Ceo model. By assuming an unre- 
alistically small phonon frequency for the intermolecular 
vibrations, we can obtain the same value of A as for the 
LE coupling. For the HI coupling, we then find that the 
resistivity our Ceo model shows saturation. We have also 
considered LE coupling in the Nbg model. The resulting 
resistivity shows a change in slope, somewhat similar to 
Fig. 1, but with a larger slope for large T than in Fig. 1. 
These results illustrates that it is possible to obtain sat- 
uration with LE coupling, but that HI coupling is more 
appropriate for describing saturation. 

Our scmiclassical treatment of the two models is closely 
related to conduction in disordered systems with diagonal 
(DD) or (ODD) off-diagonal disorder. While DD models 
can give localization, no localization was not found close 
to the center of the band in a ODD model ||2^]. This 
is consistent with the saturation seen in Fig. 1 for the 
A15 model, having ODD in the semiclassical treatment. 
For the Ceo model, the semiclassical treatment gives lo- 
calization for a large DD. The QMC treatment, however, 
takes into account that the disorder is not static but due 
to thermal fluctuations and that the scattering can be 
inelastic. In this QMC treatment wc find a lack of satu- 
ration, but we have not seen signs of localization. 

To summarize, guided by the loss of periodicity at large 
T, we have studied the effect of replacing the current ma- 
trix element by its average. Together with charge conser- 
vation, this leads to clear saturation in a model where the 
phonons couple to the hopping matrix elements. On the 
other hand, saturation does not happen in a model for 
AaCeo, where the phonons couple to the level energies. 
The issue of saturation or not saturation is only raised 
for experimentally accessible temperatures if A is large 
and Qpi is small. This is favoured by the relatively flat 



bands for Nbg {^pi small) and by the small band width 
for AaCeo (A large and flpi small). 
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